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1.  INTRODUCTION. 

Consider  a  source  of  light  and  a  linear  path  in  a  plane.  Random  objects  are  scattered 
between  the  source  of  light  and  the  path,  casting  shadows  on  portions  of  it. 

The  Shadowing  Model  was  originally  introduced  by  ChernofF  and  Dally  [1].  In  their  paper 
they  assumed  that  the  shadowing  objects  are  disks  having  centers  which  are  distributed 
in  a  rectangle  between  the  source  of  light  and  the  path,  according  to  a  two  dimensional 
Poisson  Process,  and  their  diameters  axe  represented  by  i.i.d.  random  variables.  The 
paper  applies  queueing  theoretical  techniques  to  derive  the  conditional  distributions  of  the 
lengths  of  the  shadowed  and  the  unshadowed  segments  of  the  path. 

The  shadowed  and  unshadowed  segments  of  the  path  are  analogous  to  busy  and  idle 
ptnods  in  certain  queueing  processes  with  Markov  arrivals.  Like  in  these  processes,  the 
derivation  of  the  distribution  of  the  lengths  of  unshadowed  segments  (which  are  similar  to 
idle  periods)  is  straightforward.  Chernoff’s  and  Daily’s  paper  concentrates  on  the  distri¬ 
butions  of  shadowed  segments,  and  provide  some  integral  equations  for  these  distributions. 
The  solution  of  these  equations  requires  generally  lengthy  numerical  procedures. 

Yadin  and  Zacks  [2,  3]  introduced  the  concept  of  Random  Measures  of  Visibility,  which 
are  defined  as  the  total  length  of  the  unshadowed  subsegments  of  a  specified  segment  of 
the  path,  or  the  proportion  of  the  segment  which  is  in  the  light.  In  these  papers  the 
distributions  of  Measures  of  Visibility  were  approximated  by  a  mixture  of  a  two  points 
distribution  on  {0,1}  and  a  beta  distribution  on  (0,1).  These  mixed  beta  distributions 
were  obtained  by  equating  the  first  three  moments  of  the  distributions  of  Measures  of 
Visibility  to  the  first  three  moments  of  a  mixed-beta  distribution.  In  [2]  the  path  wa s 
assumed  to  be  a  segment  of  a  circle,  centered  at  the  origin  at  which  the  source  of  light  was 
located.  The  centers  of  the  disks  were  assumed  to  be  distributed  within  an  annular  region 
in  the  circle.  The  results  were  extended  in  [3]  and  [4],  which  considered  more  general  paths 
and  fields  of  obscuring  objects.  The  methodology  was  extended  to  a  three  dimensional 
model  in  [5].  Certain  applications  were  discussed  in  [7]. 
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The  main  objective  of  the  present  paper  is  to  replace  the  mixed  beta  distributions  by 
discrete  distributions  which  approximate  the  required  distributions.  More  specifically,  we 
approximate  the  distribution  of  a  measure  of  visibility,  by  the  distribution  of  a  discrete 
random  variable,  being  the  number  of  the  proportion  of  the  unshadowed  subsegments, 
which  are  properly  located  on  some  segment  of  the  path. 

The  shadowing  model  and  the  approximated  measures  of  visibility  and  lengths  of  shad¬ 
owed  segments  are  formally  defined  in  section  2.  The  approximating  distributions  can  be 
given  in  terms  of  probabilities  of  elementary  events,  which  are  the  events  signifying  which 
of  the  N  subsegments  are  completely  unshadowed  among  the  set  of  all  N  subsegments. 
These  probabilities  can  be  obtained  from  probabilities  of  simultaneous  visibility,  namely 
the  probabilities  that  all  subsegments  in  a  specified  sets  are  unshadowed.  Simultaneous 
visibility  probabilities  were  determined  for  various  cases  in  [4],  Appendix  A  presents  the 
methodology  needed  for  the  derivations  of  probabilities  of  simultaneous  visibility. 

The  accuracy  of  the  discrete  approximation  to  the  distributions  considered  depends  on 
the  number  of  subsegments.  Assuming  a  set  of  N  subsegments,  one  has  to  compute  2N 
probabilities  of  elementary  events  from  the  same  number  of  probabilities  of  simultaneous 
visibility. 

Any  algorithm  which  is  based  on  the  probabilities  of  simultaneous  visibility  is  necessarily 
non  polynomial,  thus  is  limited  to  a  relatively  small  number  of  points,  and  may  therefore 
yield  insufficient  accuracy.  Fortunately,  the  distributions  of  the  approximated  measures  of 
visibility  can  be  obtained  by  a  polynomial  algorithm  which  is  discussed  in  section  4. 

The  results  of  a  numerical  example  are  presented  and  discussed  in  section  5.  These 
results  illustrate  the  advantages  and  limitations  of  the  algorithms  developed  in  the  present 
paper  with  respect  to  accuracy,  computer  time  and  memory  requirements. 
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2.  The  Shadowing  Model 

Let  C  be  a  star-shaped  rectifiable  curve  in  the  plane,  with  respect  to  a  point  O.  C 
represents  a  possible  path  of  a  moving  target  and  O  the  location  of  an  observer.  We  further 
assume  that  C  does  not  pass  through  O.  A  countable  number  of  objects  are  randomly 
dispersed  between  C  and  O .  These  random  objects  may  obscure  the  visibility  from  O  of 
certain  portions  on  C .  The  model  assumes  that  the  obscuring  objects  are  disks  centered 
according  to  a  two-dimensional  Poisson  point  process  of  rate  A ,  in  a  region  S ,  located 
between  O  and  C  and  neither  O  nor  C  are  intersected  by  any  such  disk.  Furthermore, 
it  is  assumed  that  the  radii  of  these  random  disks  are  i.i.d.  random  variables,  having  a 
specified  c.d.f.,  F(x) ,  concentrated  on  [0,6]. 

Let  0  be  an  interval  on  the  real  line,  representing  the  indices  s  of  points  P3  E  C  so 
that,  if  Si  <  s2,  si,s2  E  0,  then  P3l  is  on  the  left  of  P32 .  Let  Ra  designate  a  ray 
originating  at  O  and  passing  through  P3.  We  say  that  P„  is  visible  from  O  if  Ra  is  not 
intersected  by  random  disks  centered  in  S .  For  each  s  €  0 ,  let  I(s)  be  a  random  variable 
assuming  the  value  1  if  P„  is  visible  from  O ,  and  the  value  0  otherwise.  The  shadowing 
process  on  C  is  characterized  by  the  indicator  stochastic  process 

(2.1)  /=  {/(<)  :se  0}. 

The  Poisson  field  assumption  implies  that  I  is  a  semi-Markov  process  in  which,  {/(r)  : 
r  >  s}  and  {/(r)  :  r  <  s}  are  conditionally  independent,  given  {/(s)  =  l},  for  each 
3  E  0.  Let  W(s)  be  a  right-continuous  non-decreasing  function  on  0.  For  specified 
points  s',s"  G  ©  such  that  s'  <  s"  define  the  random  variable 

Bn 

(2.2)  V  =  J  I{s)dW(s). 

V  is  called  a  random  visibility  measure  on  a  segment  C ,  of  C  bounded  by  P3i  and  P3n . 
The  interpretation  of  V  depends  on  the  specific  definition  of  the  function  W .  For  example, 
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if  W(-s)  is  the  distance,  measured  along  C  between  Pa>  and  Pa,  then  V  represents  the 
total  length  of  the  portions  of  C ,  which  are  visible  from  O . 

Another  example  is  given  when  W  (s)  is  the  time  required  for  a  target  to  move  along  C 
from  Pa>  to  Pa .  In  this  case  V  is  the  total  length  of  time  in  which  the  target  is  visible 
from  O ,  while  moving  from  P a>  to  Pa» .  The  weight  of  the  shadow  to  the  right  of  P3  on 
C  is  defined  as 

(2.3)  X(s)  =W{s  +  t(s))  -W(s), 
where 

(2.4)  t(s)  =  inf  {r;  r  >  0,  I[s  +  t)  —  1}. 

Notice  that  t(s)  and  X(s)  are  random  variables.  Correspondingly,  the  shadowing  process 
on  C  is  the  stochastic  process 

(2.5)  X  =  {X(s)  :  s  E  ©}. 

The  Poisson  field  assumption  implies  that  X  is  a  Markov  process,  in  which  (-X'(r)  ;  t  <  s) 
and  {X(r);  r  >  s}  are  conditionally  independent,  given  X(s),  for  each  s  €  0.  Further¬ 
more,  all  sample  paths  of  X  are  piece-wise  nonincreasing,  having  countable  upwards  jump 
points  at  which 

(2.6)  Hm(J(s)  -  I[s  -  e))  =  -1. 

These  jump  points  are  the  left  end  limits  of  shadowed  portions  of  C.  In  the  following 
sections  we  provide  discrete  type  approximations  to  the  distribution  of  V  and  to  the 
conditional  distribution  of  X(s),  given  that  s  satisfies  (2.6). 
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3.  Discrete  Approximations  of  V  andX(s). 

Consider  the  segment  C  with  end  points  Pa>  and  Pa"  •  Without  loss  of  generality,  j 

assume  that  W(s)  is  a  c.d.f.  concentrated  on  [s',5"]. 

For  a  given  positive  integer  N ,  let  £,•  be  the  [i/N)th  fractile  of  W ,  i.e., 

(3.1)  &  =  W-1(^-)  =  inf{5  '  s>  s'  and  W(s)  >  } ,  *  =  0, 

Notice  that  £0  =  s'  and  =  s".  These  fractiles  define  a  partition  of  C  to  N  subsegments. 

The  measure  of  visibility  V  can  be  written,  in  terms  of  this  partition  as 

*  /■«* 

(3.2)  V=52 

Ui-i 

Let  P9i  (*  =  1,..  .,N)  be  the  “midpoint”  of  the  tth  subsegment,  defined  by 

(3.3)  i  =  l . N. 

Let  Jrtts(si){i  =  1, . . . ,  N)  be  a  subsegment  of  weights  S/N,  centered  at  P3i ,  0  <  <5  <  1 . 

These  subsegments  are  specified  by 

(3.4)  JN,s[si)  =  {Pa  •  5,-  ^  5  <  5,'  }  > 


where 

(3.5) 


1 


Finally,  let  lN,s{si )  be  an  indicator  random  variable,  which  assumes  the  value  1  if  JN,s(si) 
is  completely  visible. 

The  measure  of  visibility  V  can  be  approximated  by  replacing  the  function  J(s)  in  (3.2) 
with  the  indicators  Is,6{^i)-  Accordingly,  the  random  variable  V  is  approximated  by  a 
discrete  random  variable  Vjv.s  >  given  by 

vn,s  =  T,  /  IN,e(si)dW(s 

tel  •'«— 

=  Jf  52 

i-1 


6 


This  approximation  converges  uniformly  to  V  (irrespective  of  5),  with  probability  one. 
Indeed, 


(3.7) 


\Vn,c  -  v\  <  i^sup{|/v,{(s,)- 

1=1 

/(s)|;  Zi-i<s<Zi}<^K(s',s")t 


where  K(s',s,t)  is  the  number  of  shadows  on  C.  The  Poisson  assumption  implies  that 
K(s',s")  is  finite  with  probability  one,  irrespective  of  N  and  6. 

Furthermore,  for  every  realization,  Ijv,s(st)  is  a  decreasing  function  of  6.  Hence,  VN>5 
is  a  decreasing  function  of  6 ,  converging  to  Vjv.i  a.s.  as  8  — ►  1.  In  addition,  Vjv,i  <  V 
for  each  N ,  a.s..  This  means  that,  for  values  of  6  close  to  1 ,  the  convergence  of  to 
V  is  from  below. 

The  shadow-weights  X(s) ,  for  £t_i  <  s  <  are  approximated  by  Xjv,fi(st),  i  = 
1,...,JV;  which  are  discrete  random  variables  defined  as  the  proportion  of  subsegments 
JN,s(si+j) ,  on  the  right  of  which  are  shadowed,  i.e. 

(3.8) 

1  -  '-JT,  if  /(«i+m)  =  0 

for  all  m  =  0, . . . ,  N 

/ 

jf  inf{m  :  0  <  m  <  N  —  *  for  which  /(st+m)  =  l},  otherwise. 

Using  the  previous  arguments  about  the  nature  of  the  Poisson  process,  we  obtain  that 


XNts{^i)  = 


(3.9)  lim  sup  |Xjv,5(«s,)  —  X(s)|  =  0,  a.s., 

iV~°°  $<-!<«<£.■ 

for  each  6 ,  0  <  6  <  1 .  Moreover,  for  each  N  >  1 ,  and  t  =  1, . . . ,  N ,  £,_i  <  s  <  £,  , 


(3.10)  X(s)  <  Xtf'i(si),  a.s.. 

We  discuss  now  the  approximations  to  the  distributions  of  V  and  X(s),  given  by  the 
distributions  of  the  discrete  random  variables  Vjv,$  and  t  =  1, . . .  ,N}. 
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Let  H(v)  denote  the  c.d.f.  of  V.  As  shown  in  [4]  and  in  Appendix  A,  the  c.d.f. 
of  V ,  under  general  conditions  on  the  distributions  of  the  radii  of  disks,  has  two  jump 
points,  at  v  =  0  and  at  v  =  1 ,  and  is  absolutely  continuous  on  the  open  interval  (0, 1) . 
po  =  P{V  =  0},  is  the  probability  that  C  is  completely  shadowed,  and  pi  =  P{V  =  1}, 
is  the  probability  that  C  is  completely  visible.  Thus,  generally 

'0,  if  v  <  0 

(3.11)  H(v)  =  '  po  +  Jq  h(y)dy,  if  0  <  v  <  1 

<  1,  l<v 

where  h(y)  >  0,  Riemann  integrable,  and 

(3.12)  f  h{y)dy  =  1  -  pi  -  po- 

Jo 

The  c.d.f.  of  V ,  H (v) ,  is  approximated  by  the  c.d.f.  of  the  discrete  r.v.  VN>5  .  Let 

(3.13)  hpfis(k)  =  P{VNts  =  -^r},  k  =  0,l,...,N. 

The  c.d.f.  of  VN6 ,  HNiS(x ) ,  is  a  step  function,  given  by 

'  0,  if  x  <  0 

(3.14)  Hn,s{x),=  <  £y=o  hNiS{j),  if  £  <  x  <  ^  (fc  =  0, 1, . . . ,  JV  -  l) 

.1,  if  1  <  x. 

Notice  that  h^ts( 0)  is  the  discrete  approximation  to  po  and  ^^^(l)  is  an  approximation 
to  p\.  Notice  that  /^^(l)  =  p\ ,  and  that  Hn,i{ k)  >  H(jt),  for  each  N  >  1,  k  = 
0,1 , . . . ,  iV.  Moreover,  po  <  h^5(0)  for  every  0  <  <5  <  1 . 

The  distribution  of  X(s )  is  concentrated  on  the  interval  [0,  1  —  W($)].  Typically,  its 
c.d.f.  has  jump  points  at  the  two  end  points  of  the  interval,  and  is  absolutely  continuous 
in  its  interior.  We  are  interested  in  the  conditional  distribution  of  X(s) ,  given  that  P3  is 
in  a  shadow,  and  the  weight  of  the  shadow  to  the  left  of  P3  is  exactly  y,  for  some  y  >  0. 


Let  G(x;s,y)  denote  the  tail  probability  of  this  conditional  distribution.  Formally 


(3.15)  G(x;s,y)  =  P{X(s)  >  x  |  J(r )  =  0,  Vty(s)  <  r  <  s}, 
where 

(3.16)  ty(s)  =  sup{2  :  I[t)  =  1  and  IF(.s)  —W(t)  >  y}. 

The  function  G(x;  s,y)  can  also  be  expressed  in  terms  of  the  process  I  as 

G(x;s,  y)  =  P{/(u)  =  0  for  all 

(3.17)  s  <  u  <  W-1(W(s)  +  x)  j  /(r)  =  0  for  all 

ty(s)  <t  <s}. 

For  each  s  in  (t  =  1  and  y  =  jf  (v  =  0, . . .  ,t  -  1)  we  approximate 

G(x;s,y)  by 

(3.18)  GjviS(x;t»  =  ^2 1{j  :  J  -  N  x  +  *.*'). 

i 

where  gN,s{j\i,v)  is  the  conditional  probability  function 

i+3- 1  i 

(3.19)  9N,s{j] *\ v)  =  P{  (l  —  /;)/j'4-y  =  1  |  (1  —  //)  =  1}. 

l=i  1=1- V 

In  the  following  sections  we  provide  an  algorithm  for  the  computation  of  the  probability 
functions  hNiS(k)  and  gNi5(j;i,v). 
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4.  Recursive  Determination  Of  hjv,s(fc)  and  gN,s{j\  *\  *>)  • 

In  the  present  section  we  develop  a  recursive  algorithm  for  the  determination  of  the 
probability  functions  under  consideration,  for  specified  values  of  N  and  8.  For  the  sake 
of  simplicity,  we  replace  Is,s{si)  and  Jjv,$(. s,)  by  and  J;,  respectively. 

We  start  with  definitions  of  probability  functions,  which  play  a  role  in  the  algorithm. 

(i)  Fora  given  t,  *  =  1, . . .  ,N,  let  £_  (t)  [resp.  £+(t)  ]  denote  the  probability  that  the 
subsegment  J,  is  completely  unshadowed  by  disks  centered  to  the  left  [resp.  to  the 
right]  of  the  ray  R#i . 

(ii)  For  given  indices  1  <  *  <  j  <  N ,  let  £(i,j)  designate  the  probability  that  J,  and 
Jj  are  completely  unshadowed  by  disks  centered  between  the  rays  Ra;  and  Ray . 

(iii)  For  n  =  0,  —  1,  m  =  1  ,...,N  —  n  and  k  =  0,  ...,m  —  1,  let  Q(n,  m,  A:)  denote 

the  probability  that  Jn,  Jn+m,  and  exactly  k  out  of  the  m  —  1  subsegments  in 
between,  are  completely  unshadowed  by  disks  centered  between  the  rays  Raji  and 

R*n+m  ■ 

For  n  =  0,  Q(0, m, k)  is  the  probability  that  Jn,  as  well  as  exactly  k  of  the  m  —  1 
subsegments  to  its  left,  are  completely  unshadowed  by  disks  centered  to  the  left  of  RSm  . 
From  the  above  definitions  we  obtain  the  following  relationships: 

(4-1)  Q(0,l,0)  =  e-(1), 

and 

(4.2)  <3(n,l,0)  =  f(n,n  +  l),  n=  l,...,iV-  1. 

Furthermore,  the  renewal  nature  of  the  processes  implies  the  recursive  formulae: 

m—  1 

(4.3)  Q{n,m,k)  =  ^  Q{n,v,k  -  1  )Q{n  +  is,  0), 

i/=fc 

for  m  =  2, . . . , N  —  1 ,  n  =  0,...  ,N  —  m  and  k  =  1, . . . , m  —  1 .  In  addition,  since 

m—  1 

(4.4)  ^  Q(n,m,/c)  =  ^(n,n  +  m), 

k=o 
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we  obtain  that 


(4.5) 


m—  1 

Q(n, 171,0)  =  £(n,n  +  m)  —  ^  <5(n,m,A:). 

fc=i 


The  functions  £±(n)  and  £(n,  m)  for  special  cases  of  two- and  three-dimensional  shadowing 
models,  with  6  —  0,  were  previously  presented  in  [4]  and  [5].  In  Appendix  A  we  present 
the  general  methodology  of  computing  these  functions. 

The  algorithm  for  computing  the  probability  function  hn,s{k),  k  =  0, . . . ,  N ,  is  given 
by  the  following  recursive  formulae: 

Let  Pn(k) ,  for  n  =  1  and  k  =  0,1,...,  n,  denote  the  probabilities  that  k  out 

of  the  first  n  subsegments,  Ji,...,Jn,  are  completely  unshadowed.  We  compute  these 
functions  recursively,  according  to  the  formulae 


<3(0,1,0)£+(1),  *  =  1 


(4.6) 


Pi(k)  = 


1  -  Pl(  1). 

and  for  n  —  1,2, . . . ,  N , 

(  Q{0,n,n-  lK+(n), 


k  =  0 


(4.7)  Pn(k)  =  { 


k  =  n 


Pn-i[k)  +  [(5(0, n,k  —  1)  -  <5(0,n,fc)]£+(n),  k  =  l,...,n  -  1 


1  1  Pn{k), 


k  =  0. 


Finally, 


(4.8)  hff,s{k)  =  Pn{Ic),  k  =  0,  l,...,iV. 

For  the  purpose  of  determining  the  conditional  probability  function  we  intro¬ 

duce  the  function 

JQ(n,m,0),  if  n  +  m  <  N 


£+(n)  ~  £u=im  ?(n-1/)^+(n  +  if  n  +  m  =  N  +  1. 


When  n  >  1  and  n  +  m  <  N  then  the  function  q(n,m)  yields  the  probability  that  the 
subsegments  Jn  and  Jn+m  are  completely  unshadowed  by  disks  centered  between  R3  and 
R«n+m  >  while  all  the  (m  - 1)  subsegments  {  Jrt+1 Jn+m}  are  shadowed.  g(0,  m)  [resp. 
q{N  +  1  -  m,m)]  is  the  probability  that  Jn  [resp.  JN+l_m]  is  completely  unshadowed 
by  disks  centered  to  the  left  [resp.  to  the  right]  of  Ra>n  [resp.  Ra„+I_m],  while  all  the 
subsegments  to  its  left  [resp.  to  its  right]  are  shadowed. 

Finally,  the  conditional  probability  function  of  Xs,6  is  given  by 


(4.10) 


9N,s{j]i ,  u) 


q(i-  v  -  l,j  +  u  +  l)Z+(i  +  -  v  1) 

Yli=  i  +1  g(*  — 1>  —  1,/  + 1/  +  !)£+(*  +  /)£-(*  —  v  —  1) 


for  j  =  1, . . . ,  N  -  i  +  1 ,  where  £+  (N  +  1)  =  1  and  (0)  =  1 . 
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5.  A  Special  Case  Of  Linear  Path  And  Trapezoidal  Region. 

In  the  present  section  we  illustrate  the  discrete  approximations  for  a  special  geometrical 
configuration.  We  consider  a  shadowing  process  in  the  plane.  The  source  of  light  is  at  the 
origin,  O  =  (0, 0) .  The  path  curve  C  is  a  straight  line  perpendicular  to  the  ray  R0 ,  at 
distance  R  from  O .  The  line  segment  C  is  bounded  between  points  Pai  and  P a»  on  C ; 
where  Pa  is  a  point  on  the  line  C,  having  polar  coordinates  (s,iE/cos(s)) . 

The  centers  of  disks  are  randomly  dispersed,  according  to  a  two  dimensional  homoge¬ 
neous  Poisson  process  of  intensity  A,  between  two  lines,  U ,  W  ,  parallel  to  C,  at  distances 
u  and  w ,  from  O ,  respectively.  We  further  assume  that  the  radii  of  disks  are  independent 
random  variables  having  a  common  uniform  distribution  on  [0,6].  Disks  which  can  cast 
shadows  on  C  must  have  centers  within  the  trapezoid  bounded  by  U ,  W  and  the  rays  R#< 
and  Rgn ,  where 


(5.1) 


O'  =  tan  x(u  tan(s')  —  6/cos(s')), 


0"  —  tan~l{u  tan(s/;)  +  6/cos(s")). 

According  to  the  general  methodology  for  computing  visibility  probabilities  (see  Yadin 
and  Zacks  [3],  [4])  we  introduce  the  auxiliary  functions  AT+(s, t)  and  K-(s,t ),  where 
XK±(s,t)  is  the  expected  number  of  disks  centered  between  the  rays  Ra  and  Rs±t,  which 
do  not  intersect  Ra.  As  shown  in  detail  in  [7], 


(5.2)  K-(s,t)  =  K+{— s,t),  for  all  s'  <  s  <  s",  t  >  0, 


and 

(5.3) 

where 

(5.4) 


tf+^M).  *<*?(£> 5) 


K+{s,t)={  K{+\s,t),  <  t  <  »?(£,$) 


I  ^+3)(5,t),  »?(£,«)<*, 


sc 

r](x,s)  =  tan-1(tan(s)  -I - — )  -  s, 

cos(s) 
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and 


(5.5) 


K+’M  = 


sin2(f) 

cos(-s)  cos 2(s  +  t)’ 


(5.6) 


w2 .  ,  .  .  b2 

T|ton(J  +  t)-tan(*)l+g^) 


,m.3  sin;(t) 

'  b  cos2  (s  +  t) 


cos  (s  +i) 
sin(t) 


(5.7)  K<?)(3,t)  =  K?\s,rJ(-,s))  + 

tt 


w2  -  u2 


[tan(£  +  a)  —  tan(s)]  — 


w2  —  u2 


2  u  cos(s) 


Consider  the  weight  function  on  [s',  s" 


(5.8) 


VF(s)  =  T(*an(5)  —  tan(s,))> 
L 


where  L  =  i2(tan(s")  —  tan(s')) .  VF(s)  is  the  proportional  distance  along  C  of  P3  from 
Pa> ,  for  s'  <  s  <  s" .  Accordingly,  the  measure  of  visibility  under  consideration  is  the 
proportion  of  the  total  length  of  C  which  is  visible,  i.e., 


(5.9) 


V  = 


R 

L 


ds 

cos2 (s) * 


For  a  given  N  and  6 ,  0  <  6  <  1 ,  we  approximate  V  by  Vs,s ,  which  is  obtained  by 
partitioning  the  line  segment  C  to  N  equal-length  subintervals.  The  points  P3n  (n  = 
1, . . . ,  N)  are  the  midpoints  of  these  subintervals.  The  subsegments  Jn  (n  =  1, . . . ,  N)  are 
centered  around  P3n  ,  with  end-points  Pa-  and  Pa+  of  distance  6/2N  from  P3n  . 

The  functions  £±(n)  and  £(i/, n),  0  <  u  <  n  <  N ,  defined  in  the  previous  sections  are 
given,  in  terms  of  the  K -functions,  by  the  following  formulae: 


(5.10) 


{-(n)  =  exp{-A(r(0',s„)  -  K. (»„,»„  -  8')]}. 


(5.11) 


{+(»)  =  exp{-A[T(s„,fl")  -  K+(siJ"  -  4)]}, 


and 


(5.12) 


{(k,»)  =  exp{-A[r(«„,s„)  -  lf+(i+, 


where 


10^  —  ti^ 

(5.13)  T{0u02)  =  - - - (tan(02)  -  tan(0i)), 

is  the  area  of  the  trapezoid  bounded  by  U ,  W ,  R$t  and  R$a . 


6.  A  Numerical  Example. 

As  in  the  previous  section,  we  consider  a  linear  curve  with  a  trapezoidal  region  for  disk 
centers.  The  radii  disks  are  independent  random  variables  having  a  uniform  distribution 
on  the  interval  (0,0.4) .  The  path  C  is  a  horizontal  straight  line  of  distance  R  =  1.0  from 
the  origin.  The  horizontal  boundaries  of  the  strip  S  are  at  distances  u  =  .4  and  w  =  .6, 
from  O . 

The  segment  £  is  the  interval  between  s'  =  —  1  and  s"  =  1  on  C.  The  intensity  of  the 
Poisson  field  is  A  [disk  centers/unit  squared  of  area].  In  Table  6.1  we  present  the  c.d.f. 
of  Vj v,s,  for  N  =  40,50,  6  =  0,1.0  and  A  —  10,20.  The  c.d.f.  Hn,6{x)  is  tabulated  for 
x  =  0(0.1)  1~.  #*i5(0)  is  the  discrete  approximation  to  the  value  of  p0.  1  —  Hn,i{  1— ) 
is  the  discrete  approximation  to  the  value  of  p\ .  In  addition,  we  provide  in  Table  6.1  the 
mixed-beta  approximation  to  the  c.d.f.  # ( x ) .  The  mixed-beta  c.d.f.  is  the  function 

'  Po,  x  =  0 

(6.1)  H*(x)  -  <  po  +  (1  -  po  -  px)Ix{a,p),  0  <  x  <  1 

.  1,  1  <  x 

where  Ix(a,/3)  is  the  incomplete  beta-function  ratio.  The  values  of  po,  a  and  (3  are 
determined  by  equating  the  first  three  moments  of  V ,  which  can  be  numerically  determined 
(see  Yadin  and  Zacks  [3]),  to  the  first  three  moments  of  H*(x).  As  seen  in  Table  6.1, 
-flso.i (z)  <  #40,1  (a:),  for  all  0  <  x  <  1.  Indeed,  Hn,i(x)  converges  monotonically  to 
#(x)  from  above.  Similarly,  #50,0  (0)  <  #40,0  (0)  in  all  the  two  cases  of  A  =  10  and  20. 
Indeed,  #n,o( 0)  >  po  for  all  N ,  and  #v, o(0)  — ^ ►  po  as  N  — *  00 .  In  Table  6.1  we  see  that, 
for  each  '  the  values  of  the  mixed-beta  c.d.f.  H*(x)  satisfy 

#50,0(1)  <  #*(x)  <  #5o,i (x) ,  x  =  .2(.l).8. 

This  indicates  that  the  mixed-beta  approximation  #*(x)  is  apparently  close  to  the  true 
c.d.f.  #(x),  over  a  wide  range  of  x.  In  Table  6.2  we  compare  the  moments  of  #jv, o(^)> 
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Hn.x  (z)  and  S' *  (x) .  We  see  that  the  moments  of  Hn,o{x)  are  very  close  to  those  of  S*(x) . 
Recalling  that  the  first  three  moments  of  H*(x )  are,  by  definition,  the  corresponding 
moments  of  V ,  we  conclude  that  the  moments  of  Hn) o(x)  are  close  to  the  true  ones. 
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A 

X 

N  =  40 

N  =  50 

N  =  40 

o 

to 

II 

H’lx) 

0.0 

0.00081 

0.00072 

0.00225 

0.00167 

0.0000 

10 

0.1 

0.00461 

0.00501 

0.01032 

0.00962 

0.0021 

0.2 

0.02325 

0.02461 

0.04286 

0.04036 

0.0305 

0.3 

0.07009 

0.07315 

0.11159 

0.10657 

0.0942 

0.4 

0.15805 

0.16336 

0.22352 

0.21624 

0.1967 

0.5 

0.29110 

0.29864 

0.37336 

0.36519 

0.3339 

0.6 

0.45836 

0.46724 

0.54266 

0.53547 

0.4954 

0.7 

0.63657 

0.64548 

0.70731 

0.70265 

0.6645 

0.8 

0.79584 

0.80301 

0.84177 

0.83984 

0.8190 

0.9 

0.90577 

0.91018 

0.92848 

0.92824 

0.9325 

Pi 

0.02620 

0.02520 

0.02140 

0.02140 

0.0214 

0.0 

0.01588 

0.01448 

0.03489 

0.02796 

0.0084 

20 

0.1 

0.06019 

0.06440 

0.10746 

0.10287 

0.0787 

0.2 

0.18754 

0.19516 

0.27790 

0.26839 

0.2372 

0.3 

0.36854 

0.37822 

0.47959 

0.46829 

0.4278 

0.4 

0.56755 

0.57417 

0.66784 

0.65805 

0.6121 

0.5 

0.73678 

0.74463 

0.81370 

0.80708 

0.7668 

0.6 

0.86270 

0...8802 

0.90931 

0.90579 

0.8803 

0.7 

0.93964 

0.94263 

0.96245 

0.96102 

0.9515 

0.8 

0.97839 

0.97974 

0.98712 

0.98672 

0.9869 

0.9 

0.99400 

0.99449 

0.99651 

0.99646 

0.9983 

Pi 

0.00069 

0.00054 

0.00046 

0.00046 

0.0005 

,  =  10 


A  =  20 


6  = 

=  0 

o 

II 

N  =  50 

Kfml 

0.5933 

0.3950 

0.2843 

0.2842 

0.2167 

0.2165 

0.1726 

0.1725 

0.1425 

0.1423 

0.1209 

0.1207 

0.1051 

0.1048 

0.0930 

0.0928 

0.0837 

0.0834 

0.3556 

0.1031 

0.3556 

0.1630 

0.0868 

0.0867 

0.0513 

0.0512 

0.0329 

0.0329 

0.0225 

0.0225 

0.0162 

0.0162 

0.0122 

0.0122 

0.0095 

0.0095 

0.0077 

0.0076 

N  =  40 


0.1402 

0.1140 

0.0957 

0.0825 

0.0726 

0.0650 


.3043 
.1278 
.0638 
0.0360 
.0222 
.0148 
0.0104 
0.0077 
0.0059 
0.0047 


N  =  50 


0.5574 

0.3565 

0.2493 

0.1860 

0.1458 

0.1188 

0.0999 

0.0861 

0.0758 

0.0679 


0.3139 

0.1340 

0.0677 

0.0385 

0.0239 

0.0160 

0.0113 

0.0084 

0.0064 

0.0051 


H*(x) 


0.5933 

0.3949 

0.2840 

0.2158 

0.1710 

0.1400 

0.1177 

0.1011 

0.0884 

0.0785 


.3556 
.1629 
0.0865 
0.0508 
0.0321 
0.0215 
0.0151 
0.0109 
0.0082 
0.0063 


We  conclude  this  section  with  a  numerical  example  of  the  conditional  distribution  of  the 
remaining  shadow  length,  given  that  v  subsegments  to  the  left  of  the  i  th  subsegment  are 
in  shadow.  In  Tables  6.3  and  6.4  we  present  the  conditional  p.d.f.  (4.10),  for  the  case  of  a 
trapezoidal  region,  with  R  =  1,  u  =  .4,  w  =  .6,  s'  =  -1,  s"  =  1,  6  =  .2  and  \  =  25,50. 
The  interval  [s'js"]  is  partitioned  into  N  =  20  subintervals.  Table  6.3  presents  these 
distributions  for  5  =  0  and  Table  6.4  for  5  =  1.  For  the  sake  of  space  saving,  we  present 
these  conditional  PDF’s  for  t  =  10;  j  =  1,2, ...,11  and  u  =  0,1,2  and  3  only.  Since 
N  =  20,  j  =  11  corresponds  to  the  case  that  all  the  10  subsegments  to  the  right  of  J\0 
are  shadowed.  We  see  that  when  5  =  0,  the  modes  of  all  the  conditional  distributions  are 
at  j  =  1;  which  means  that  it  is  most  likely  that  the  eleventh  subsegment,  Jn ,  is  in  the 
light.  This  is  not  necessarily  the  case  when  5  =  1. 


Table  6.3  Values  of  the  Conditional  P.D.F.  gij'ji,^),  for 
i  =  10,  in  a  discrete  approximation ,  with  N  =  20,  5  =  0. 
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10 
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10 
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■ 
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11 
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9 
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.2817 

.1904 
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Appendix  A.  Visibility  Probabilities  For  Poisson  Fields. 

We  present  here  the  basic  methodology  of  determining  visibility  probabilities,  under 
Poisson  shadowing  processes,  leading  to  general  formulae  of  f±(n)  and  £(n,  m) .  For  more 
details  and  discussion  see  Yadin  and  Zacks  [2],  [3],  [4]. 

We  assume  that  a  countable  number  of  disks  are  randomly  dispersed  in  a  region  in  the 
plane  according  to  a  Poisson  random  field  model.  More  specifically,  let  ( 9,p )  be  the  polar 
coordinates  of  the  location  of  a  center  of  a  disk,  with  respect  to  an  origin  O ,  and  let  r  be 
the  radius  of  a  disk.  Let 

(A.l)  $  =  {(0)P,r)  :  (9,p,r)  G  5}, 

be  the  sample  space  of  possible  disks,  and  let  T  be  the  Borel  <7-field  of  subsets  of  5  .  Let 
T  G  T  and  N{T}  be  the  number  of  disks  having  coordinates  in  T.  The  Poisson  field 
model  assumes  that,  for  any  t  >  1  and  any  partition  {Ti, . . . ,  Tt}  of  S  ,  iV{T’1}, . . . ,  N{Tt} 
are  independent  random  variables,  having  Poisson  distributions  with  means 

(A. 2)  v{Ti}  =  A  HI  h[Q ,  p)dF{r  \  O,p)d0dp,  *  =  1 . t, 

Ti 

where  0  <  A  <  oo  is  the  intensity  of  the  field;  F(r  |  0,  p)  a  conditional  c.d.f.  of  r,  given 
( 0,p )  and  h(6,p)  is  a  bivariate  density  function  on  (— n,  7r]  x  (0,oo). 

Let  P  =  (s,p)  be  a  point  in  the  plane,  and  let  R(s)  be  a  ray  from  O,  passing  through 
P.  Let  C/_(s)  [resp.  C/+(s)]  be  all  points  in  S ,  on  the  left  [resp.  right]  of  R(s),  which 
intersect  R(s),  i.e.,  if  d(0,p;R(s))  denotes  the  distance  between  P  =  ( 0,p )  and  R(s)  then, 

(A.3)  (7_(s)  =  {(0,p,r)  :  (0,p,r)  G  5,  9  <  s,  and  d(0,p;R(s))  <  r}. 

Similarly, 

(A.4)  J7+(s)  =  {(0, p, r)  :  (0,p,r)  G  S,  9  >  s,  and  d(0,p;R(s))  <  r}. 
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Let  A M+(s)  [resp.  AM_(s)]  be  the  expected  member  of  disk  on  the  right  [resp.  left]  of 
R(s)  which  intersect  it,  i.e., 

(A. 5)  M±(s)  =  JJJ  h{6,p)dF{r  |  0,p)]d9dp. 

U±(a) 

The  probability  that  R(s)  is  not  intersected  by  disks  on  its  left  [resp.  right]  is  £_(s)  [resp. 
!+(«)].  where 

(A.6)  |±(s)  =  exp{-AM±(s)}. 

The  visibility  probability  of  P  =  (s,p)  is  £(s)  =  £+(.s)£_(s). 

If  C  is  a  segment  of  a  curve  in  the  plane,  and  Pn  =  (sn,pn),  n  =  are  “mid¬ 

points”  in  a  partition  of  C  into  N  subsegments,  then,  in  the  notation  of  the  previous 
sections  £±(n)  =  £±(sn) ,  n  =  1 .,N . 

Let  Pi  =  (s i , pi )  and  P 2  =  (32,^2)  are  two  distinct  points  and  R(s,),  1  =  1,2  are  the 
corresponding  rays  from  the  O  through  the  points.  Let  U(si,  s2)  be  the  set  of  all  possible 
disks  in  S ,  between  R(si)  and  R(s2)  which  intersect  either  R(si)  or  R(s2),  i.e., 

(A. 7)  U(si,s2 )  =  {{9,P,r)  :  ( 9,p,r )  e  S  and  d(0,  p;  R(si))  <  r  or  d(0,  p\  R(s2))  <  r}. 

The  expected  number  of  disks  between  R(si)  and  R(s2)  which  intersect  either  one  of  these 
rays  is  XM(si,  s2) ,  where 

(A. 8)  M(si,s2)  =  Iff  h(9,p)dF{r\e,p)dedp. 

u{*  i,*3) 

Finally,  the  probability  that  neither  R(si)  nor  R(s2)  are  intersected  by  disks  centered 
between  them  is  ((si,s2)  =  exp{— AM(si,32)}-  The  value  of  the  function  £(711,712)  for 
points  Pi  and  P2  on  C  is  given  by  £(7i!,n2)  =  £(-Si,32). 
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